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Abstract 

The semiclassical path integral (SPI) method has been applied for studying spin relaxation in a 
narrow 2D strip with the Rashba spin-orbit interaction. Our numerical calculations show a good 
agreement with the experimental data, although some features of experimental results are not clear 
yet. We also calculated the relaxation of a uniform spin-density distribution in the ballistic regime 
of very narrow wires. With the decreasing wire width, the spin polarization exhibits a transition 
from the exponential decay to the oscillatory Bessel-like relaxation. The SPI method has been also 
employed to calculate the relaxation of the particularly long-lived helix mode. A good agreement 
has been found with calculations based on the diffusion theory. 



PACS numbers: 72.25.Rb, 71.70.Ej, 72.25.Dc, 03.65.Sq 



I INTRODUCTION 



The spin relaxation rate is an important spin transport parameter. Recent calculations 
and measurements of this parameter in semiconductor systems have to a great extent been 
motivated by numerous ideas of spintronic applications. [1] In view of these applications, 
as well as from the fundamental point of view, one of the most interesting problem is the 
spin relaxation in quantum dots (QD) and quantum wires (QW). In zincblende semicon- 
ductors at low temperatures the spin lifetime is mainly determined by the Dyakonov-Perel 
m mechanism associated with spin-orbit effects. In systems with restricted dimensions this 
relaxation mechanism is strongly suppressed, as has been calculated in the case of QD [3|] 
and QW The physics of such a suppression in QW became clear from the an- 

alytical solution of the diffusion equation for nonuniform spin distributions confined in a 
wire.j4| Surprisingly the suppression starts when the width w of a wire becomes less than 
the characteristic length L^o of the spin orbit interaction. In typical zincblende semiconduc- 
tor systems it varies from several thousands A up to several microns and can be much larger 
than the electron mean free path I. Hence, such a spin lifetime enhancement can not be 
considered as a manifestation of the motional narrowing effect when a restricted geometry 
of the system imposes the upper limit on the mean free path. Indeed, recent measurements 
j3] have demonstrated that the spin lifetime starts to increase already at ^ 10 /. On 
the other hand, the observed slowdown of the spin relaxation appears to be not so strong, 
as expected from the theory. To understand such a behavior, one has to take into account 
that in experiments j?! the measured parameter is the relaxation time of a particular spatial 
spin distribution, rather than of an individual electron spin. At the same time, as shown 
in Ref. U only two kinds of spin distributions have very long lifetimes in narrow 2D wires. 
The first one corresponds to a polarization which is homogeneous along the wire with spins 
oriented in the plane of a 2D electron gas (2DEG) and perpendicular to the wire axis. The 
second distribution is a nonuniform helix mode with the wavelength determined by Lgo- As 
it will be pointed out below, none of these distributions have been excited by an incident 

Q 

light beam in the experiment [7|]. 

In order to interpret experimental data we will analyze relaxation of various spin distri- 
butions. We will study diffusive, as well as ballistic regimes of electron motion in the wire. 
The path integral method previously applied to QD [3] will be employed to calculate the 
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spin relaxation in a wide parameter range, including the ballistic regime w ^ I, and at time 
intervals less than the electron momentum relaxation time. 

The article is organized by the following way: sec. II gives an introduction to the path 
integral method; in sec. Ill the spin relaxation of a homogeneous spin distribution is calcu- 
lated and a comparison with the experiment is given; in sec. IV some analytical results are 
presented useful for understanding the spin relaxation behavior in the ballistic range; sec. V 
is devoted to an analysis of long-lived helix spin distributions, beyond the diffusion theory 
of Ref. y. The conclusion is presented in sec. VI. 

II SEMICLASSICAL PATH INTEGRAL APPROACH 

The semiclassical path integral (SPI) formalism has been used to study the spin relaxation 
and spin current transmission in systems with the Rashba spin-orbit interaction (SOI) 

n Q 

18|, 19| . The Hamiltonian of these systems can be divided into two parts 

H = Ho + Hr, (1) 

where Hq contains the kinetic and potential energies of an electron in a 2DEG. The second 
part Hji = a{p X (t) ■ z represents the SOI, where a is the spin-orbit coupling constant, 
p denotes the electron momentum, cr stand for Pauli matrices, and z is the unit vector 
perpendicular to the 2D sample. Within the SPI method the spin-orbit interaction Hn gives 
rise to spin precession of a particle moving along a classical trajectory. A characteristic 
length determining this precession is given by Lgo = where m* is the effective mass of 
the particle. In real semiconductor materials, the energy ratio Hr/Hq can reach 1/10, like 
in the InSb sample. [lo| But even for such a ratio, Hn is still small compared with Hq. In 
such systems whose characteristic size, or the electron mean free path, are smaller than L^o, 
the electron dynamics is not affected strongly by its spin dynamics, so that in the leading 
approximation classical trajectories are determined by Hq. 

For a free electron moving along a straight trajectory 7 of length /, the dynamics of its 
spin state is governed by the evolution operator U in the path integral formalism 



U = exp 



Hnit) dt 



7 



exp 



—I- — b ■ cr 



(2) 



where b = z x p/\p\. This operator represents simply the spin rotation. Note that the 
semiclassical approximation has been employed in Eq. ([2]), since only the classical paths 



were taken into account in the integral, This approximation is vahd for systems whose 
size is much larger than the de Broglie wavelength of electrons. The term 'semiclassical' 
here is referred to the semiclassical (saddle point) approximation taken in the spin evolution 
operator in Eq. (2). Though this terminology is used, the following analysis within the 
SPI method is classical one. This method is equivalent to the classical Boltzmann equation 
modified to take into account spin dynamics. Such an equation was used already in Ref. ^ 
and many times since. 

If an electron collides with impurities or boundaries — 1 times, its trajectory 7 will 
consist of straight segments 

7 = 7n^ + ■ ■ ■ + 7i + ■ ■ ■ + 72 + 7i, 
The corresponding spin evolution operator becomes a product 

U, = U„^^---U,^---U,,U,,, (3) 

where the individual operators 

U^^ = exp 



L 



-'so 



= 1 COS (^) - tibj ■ a) sin ( ^) (4) 

along different straight segments do not commute with each other. Each trajectory 7 at the 
time t is uniquely determined by the particle initial coordinate and momentum. Since in 
the semiclassical approximation the particle momentum is equal to the Fermi momentum 
PFi the trajectory will depend on its initial angle, while its total length is simply tvp- At 
the time t the spin of a particle moving along the trajectory 7 is determined by s'^^\t^ ctq) = 
(f/~^crf/^)o.g, where angular brackets denote averaging over the initial spin state ctq. Using 
this expression one may calculate the evolution of the spin state in systems of various 
geometries, with or without elastic impurity scatterers, as shown by several examples in 
Ref. |3|. In order to determine the evolution of a given particle distribution, one must 
perform a statistical average of the above expression over initial trajectory coordinates and 
angles, which we will denote by the trajectory label 7. Below, we will assume that particle 
initial positions and momentum directions are uniformly distributed. In this case, given a 
small area D of 2DEG, the z projection of the spin polarization Pz{t) within this area at 
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time t is determined by the average 

n(t) = ^$^3W(t,ao), (5) 

where £>) = ^^(^i) is the number of trajectories reaching the area D at time t, with 
sP'\t, (To) denoting the z-component of the electron spin. According to the above definition, 
st^\t, (To) varies within [—1, 1]. Hence, the maximum value of |-P2(t)| is 1, which corresponds 
to all electrons in D being aligned in z-direction. 

To apply the SPI method numerically, a large number of electrons are initially randomly 
distributed in the channel with uniform or helix spin configurations as explained below. 
Each electron moves straight before collision with impurities. The distance between two 
collisions follows the well known exponential distribution of free paths. In the following, we 
assume that the channels have smooth boundaries on which the electron reflection is specular. 
In the diffusive regime (as in the experimental sample in Ref. [tI), the spin relaxation 
behavior under this assumption is the same as in the case of non-smooth boundaries, because, 
even when electron trajectories are not randomized by the smooth boundary, they will be 
immediately randomized by the impurities near the boundaries. In the ballistic regime, the 
relaxation behaviors in systems with smooth and non-smooth boundaries are different. Here 
we focus on the simple example of specular reflection. Once the boundary roughness of a 
ballistic sample is known, the extension to the non-specular case is straightforward. 

Ill RELAXATION OF UNIFORM SPIN MODES 

The spin relaxation times obtained in the experiments of Ref. [3] were measured in a 2D 
n-InGaAs channel of the length L = 200 fim and the width w = 0.42 ~ 20 um. The SOI in 
the sample is dominated by the Rashba coupling. In the notation of Ref. [7| it corresponds 
to the characteristic length /^p ~ 1 /xm, which is related to the above defined spin rotation 
length Lso as L^o = '^hp- 

The sample is characterized by the electron mean free path / = 0.28 /im, the momentum 
scattering time tm = 0.76 ps, and its Fermi velocity can accordingly be estimated a.s vp = 
0.28/im^0.76 ps ~ 0.37/im/ps. For the carrier concentrations Ug = 5.4 ~ 7.0 x 10^^ cm~^ 
used in Ref. 7(], the de Broglie wavelength Xf = ^/I'k jug of electrons in the 2DEG is 
around 30 ~ 34 nm. The sample was patterned along various crystallographic directions 



and electron spins have been optically oriented parallel to the growth direction [0, 0, 1]. The 
relaxation times measured in Ref. [3] are replotted by the circles and the stars connected 
by the blue and the green curves in Fig. 1. 

Since the width range 0.4 fim < w < 20 fim used in our calculations is much larger than 
the de Broglie wavelength A/, the quantum effects are negligible and the validity of the SPI 
approach is justified. With the above experimental parameters, the SPI calculations are 
represented in Fig. 1, where the inset shows the relaxation curves Pz{t) for three channels 
of different widths. All electron spins were initially aligned in 2;-direction. The relaxation 
time Ts can be determined by a fitting of these -P^(t) curves with the exponential function 

P,{t) = Aexp{-t/T,) +c. (6) 

For example, the (red) solid curve in Fig. 1 represents the relaxation time of 1.2 x 10^ 
electrons in channels of different widths w's. A comparison with the experimental data 
(circles and stars) leads to following conclusions: 

(i) At large widths {w > 15 /im), the electron spin can be regarded as relaxing in bulk 
systems. In the experiments in Ref. |7|, Isp was estimated to be 1.0 ± 0.1 fim, corre- 
sponding to Lso = 2.0 ±0.2 yum. This experimental uncertainty results in = 9.7±2.1 
ps, when calculated by the SPI method. However, each obtained from the SPI 
method agrees very well with that determined by the analytical expression of DP re- 
laxation Tg = Ll^/{Avpl) for boundless systems. Thus, if the experimental samples 
are governed by pure Rashba Hamiltonian, as in our calculation, these samples most 
likely have Lso = 2.19 /im. This value is used in our SPI simulations to obtain the red 
and black curves in Fig. 1. 

(ii) For intermediate widths (1.4 /im< w < 15 /im), there is no an analytical expression 
for Ts to compare with. The SPI result deviates slightly from the experimentally 
measured r^. The maximum deviation is around 3 ps for [1, 1, 0] sample and 4 ps for 
[1, 0, 0] sample at ti? = 5 yum. The calculated is closer to the of the [1, 1, 0] sample. 

(iii) For small widths {w < 1.4 /im), the experimentally measured saturates at 28 ps 
for [0, 0, 1] sample and 20 ps for [0, 1, 1] sample. It is assumed in Ref. [7| that this 
saturation might be related to other mechanisms, like the bulk inversion asymmetry. 
However, the calculated in Fig. 1 is also bounded by a maximum value around 24 



ps, although in our calculation only the Rashba Hamiltonian was considered, without 
any additional mechanisms involved. 

An important factor affecting the interpretation of the experimental data is how to deter- 
mine the relaxation time from the function Pz{t)- The solid curve in Fig. 2 is an example 
of Pz{t) in a channel with w = 0.1 fim and / = 0.3 /im. At first sight it looks like an expo- 
nential function to be fitted with a relaxation time in Eq. (j6]). But a closer look shows 
that it is not a pure exponential function. Indeed, if we gradually increase / by reducing 
the number of impurities in the channel, the monotonically decreasing Pz{t) in Fig. 2 will 
transform to an oscillatory function. In the extreme case of an infinitely thin impurity-free 
channel, we shall prove in the next section that the evolution of Pz{t) will follow the Bessel 
function 



This analytical formula is depicted by the smooth (red) dashed curve in Fig. 3. A corre- 
sponding result of a numerical SPI simulation is plotted as a (green) rugged solid curve in 
the same figure. 

During transition from the diffusive to the ballistic regime, Pz(t) will undergo a crossover 
from an exponential to a Bessel function. In principle, it is meaningless to use an exponential 
function to extract from such a crossover function, especially when it is far from an 
exponential behavior. But if one would like to carry out this procedure, the so obtained 
will depend on the choice of the parameters A and c in Eq. ([6]): 

(A) If A = 1 is chosen, Eq. ([6]) can precisely fit the real initial polarization Pz{t) = 1 at 
t = (red dotted curve in Fig. 2). 

If A ^ 1, Eq. ([6]) can provide a better fitting to Pz{t) in a wider range of times at 
t > (blue dashed curve in Fig. 2). On this reason, such a choice of A seems to be 
more appropriate. 

(B) Further, if c 7^ the fitted values of c and will be strongly dependent on the 
observation time cutoff. The reason is that usually the tail of Pz{t) is oscillating, if 
the system within the considered range of times is not in the diffusive regime. The 
closer the system to the ballistic regime, the larger is the oscillation amplitude. The 
Bessel function in Eq. (I7l) for 'pure' ballistic regime has the largest amplitude. If the 
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non-oscillating Eq. is used to fit an oscillating Eq. ([7]) truncated at some cutofi^, 
the fitted Ts and c will depend on the cutoff. The corresponding uncertainty of will 
decrease with an increasing observation time. 

In the experiments the width of the channel varies between w ~ 1.5/ and 70/. Since 
at the smallest w the system is not far from the ballistic regime, the difference between 
Pz(t) and the exponential function should be observable. Indeed, the value of fitted by 
Eq. with A = 1 and c = (black dash-dotted curve in Fig. 1) is somewhat distinct 
from Ts a.t A ^ 1 and c 7^ (red solid curve in Fig. 1). Since our observation time is 
sufficiently long, the fitted value of c is close to zero. A disagreement produced by different 
fitting procedures will become more remarkable when the system approaches the ballistic 
regime with strongly oscillating Pz{t)- Hence, when comparing Ts's obtained by different 
research groups, it is important to know the whole set of the fitting parameters {A, c, and 
observation time). Even when the same Pz{t) curve is considered, the reported r^'s could be 
different. One more problem with the fitting procedure is that even in the diffusive regime 
the evolution of the spin polarization not necessarily follows the exponential behavior with 
a single relaxation time. For example, a homogeneous distribution is not an eigenstate 
of the diffusion equation in a 2D channel. Therefore, as shown in Ref. IgI, edge states can 
contribute to the Pz{t) evolution with the relaxation time different from that of the bulk 
eigenstate. The weight of edge states increases with decreasing w. 

In regime (ii), the experimental data deviate slightly from the SPI calculations with 
a maximum difference ^ 3 ~ 4 ps at w ~ 5 fim. This discrepancy is too large to 
be attributed to different fitting procedures. One of the explanations for such a behavior 
might be a specific role of long-lived edge states. The lifetime of such modes depends on the 
boundary conditions [6]. Our SPI calculations assumed a specular refiection of electrons from 
hard wall boundaries of the wire. Probably, the experimental situation in Refiri corresponds 
to other boundary conditions which give rise to the edge states with larger r^. This problem 
requires a more thorough analysis. 

In regime (iii), the relaxation time goes to a finite value at w — in both experimental 
and SPI calculated plots in Fig. 1. For a homogeneous spin distribution along the channel, 
the diffusion theory [4] also predicts a saturation of at ^ 0. The saturated value should 
be twice of the bulk DP spin relaxation time. With the experimental bulk value rs=11.5 ps, 
one expects rs=22.8 ps at w = 0. Experimental and SPI curves at Fig.l are not far from this 
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value, although the diffusion approximation fails at ~ Z. At the same time, one should not 
forget that in a narrow channel the time evolution of the spin polarization strongly deviates 
from the exponential function. On this reason, in regime (iii) Tg can not be a representative 
parameter to describe the spin relaxation. 

IV BESSEL RELAXATIONS IN BALLISTIC CHANNELS 

Depending on the ratio between the channel width and the electron wavelength, one 
encounters two limiting cases. If the wire carries only one propagating channel, we have 
effectively a ID situation. In the opposite limit, if the width of the wire is much larger than 
the electron wavelength, semiclassical electrons are able to move in both x- and t/-directions. 

Therefore, the system is two-dimensional, even if geometrically the channel is narrow, with 
w much less than other characteristic lengths, such as Lgo and /. Below, we will consider 
the evolution of electron polarization in the ballistic regime for these two limiting cases. 

ID ballistic channels 

Let us consider a ID impurity-free channel where at i = spins of all electrons are aligned 
in ^-direction. Since impurities and the electron-electron interaction are absent, electrons 
can only move in +x or —x directions along the channel axis with a constant velocity. The 
spins of all these electrons will rotate simultaneously along different geodesies connecting the 
north and south poles on the spin sphere. A 27r spin rotation takes place when an electron 
passes a distance ttL^q during a time period -kLso/vf- Therefore, the angular frequency of 
this rotation is IvpjLsoi the same for all spins. The spin polarization at any place in the 
channel will then evolve according to 



and oscillate without any amplitude decay, as shown by the blue dash-dotted curve in Fig. 




(8) 



3. 
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2D ballistic channels 



Now we suppose that the channel is a 2D thin balhstic wire where all electron spins 
are initially aligned in the z-direction, as in the ID case. Given an observation point, say 
Po{xo, Uo) in Fig. 4, the polarization Pz{t) at Xq at time t is the average of the spins of all 
electrons which will arrive at this moment at the Xq cross section. These electrons can arrive 
through a straight trajectory pTpo of length I = vpt, or through different zigzag trajectories 
of the same length, as the path P2P0 in Fig. 4. 

As follows from Eq. Q, the spin state of an electron running along the zigzag trajectory 
P2P0 in the inset of Fig. 4 will evolve according to the spin evolution operator 



P2V0 



exp —I 



Ir 



.. exp —I 



exp —I 



^ .Zn-cr , 

P2V0 ~'~ I ^2 



1 ^2'7- , 

1 - t— h. 



1-i 



licr 

Lso 



(9) 



Making the above expansion up to the linear in w term, the operator can be written 



as 



P2P0 



exp —t 



I ■ a- 



Ls 



1 - i 



I 



with the vector I pointing from p2 to po- Equation ([9]) indicates that the spin evolution of 
an electron moving along the zigzag trajectory p^Po is approximately the same as that of an 
electron drifting along the shorter straight line with a drift velocity vpx/l slower than 
vp, where / = |Zi| + IZ2I + ••• + |Z„| and x = \li + 12 + ... + ln\- As discussed in the previous 
subsection, if an electron moves a distance x, its spin will rotate the angle IxjLgo in the spin 
space. If initially this spin is aligned along the z-direction, its z-component will become 



s.(x) = cos (|^]+ 0(^1^). 



(10) 



To determine how many electrons will contribute to Pz{t), let us uniformly divide 
the channel axis into small intervals [^j, of length e separated by points C,i with 
i = 0, 1, 2, .... Let 6{^i) be the outgoing angle of an electron at ^j. This angle is cho- 
sen so that when the electron travel a zigzag path of the length / = vpt, its drift length will 
be X. Hence, this electron will arrive at po at time t. If outgoing angles are isotropically 
distributed, the number of such electrons within [^i, is proportional to the spanned 
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angle W^(Ci) = — For a given small interval e = ^j+i — this angle is related 

to X by 



W{x) = 9{x) — 9{x + e) = arccos j — arccos ^ 



2 



where the outgoing angle 6{x) of the trajectory along p2Po is the same as the angle of P2P3 
(see Fig. 4). According to Eq. the spin polarization Pz{t) at pq will be contributed from 
electrons traveling from different initial locations x: 



p W{x) Sz{x) dx 

m = , (12) 

pW{x) dx 







where p is the line density of electrons along the channel axis. Inserting Eq. (fTOj) and (fTTj) 
into Eq. ([12]) yields 



cos - — dx 



-'so 



dx 



Jo ( ]+0' 



Lso ) V-^so 

which at the w ^ limit is the Bessel function in Eq. ([7]). Hence, uniformly polarized spins 
in a ballistic narrow channel will relax to the zero polarization through a Bessel function. 
This phenomenon is in contrast to our conventional intuition that a relaxation is a mono- 
tonically exponential process. It is worthwhile to note that the Bessel-like spin dynamics 



also takes place in other SOI systems. 



We conclude that the spin relaxation dynamics even in a very thin 2D channel is re- 
markably different from that in a ID channel, as can be seen from a comparison of Eq. ([7]) 
with Eq. ([8]). It can be understood from the fact that no matter how narrow the width 
of a 2D channel is, it contains a large number of electrons moving along various zigzag 
trajectories bouncing between two channel boundaries. These trajectories give a significant 
contribution and change the dynamics of -Pz(t) from a sinusoidal oscillation in ID systems 
to a Bessel-function decay in 2D systems. 
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V RELAXATION OF HELIX SPIN MODES 



In previous sections, all initial electron spins were polarized along the 2;-axis. The re- 
laxation dynamics of such spin configuration does not change dramatically at small w (the 
maximum only reaches 28 ps in Fig. 1.), in agreement with the experiment, j?!. In fact, 
this behavior is expected from the analysis of eigenstates of the spin diffusion equation. As 
shown in Ref. 14, the only homogeneous eigenstate is \I^o iii Fig- 5(a), which has all spins 
polarized in the ^/-direction and whose relaxation time strongly increases in narrow wires. 
In addition to this mode there are two nonuniform slowly relaxing eigenstates. 

These two long-lived eigenmodes exist in a 2D channel where w -C Isp- Using a per- 
turbation method with respect to w/lsp, one can solve a diffusion equation and obtain its 
unperturbed eigensolution 

'4'M,k,m{x, y) = exp{ikx)xm{y)'^M, (13) 

with the eigenvalue (relaxation rate) 

^M,k,m = Di-nm/df + D{k - Ml~],f. (14) 

Therein, \E'a/ are the eigensolutions of the momentum operator Jy with eigenvalues M = 
0,±1 and X2n{y) = cos{27iyn/d), as well as X2n+i(y) = — sin [???/( 2n + l)/d\. In Fig. 5(a), 
will construct the helix eigenmodes with the wave vectors k = ^1/lsp- We note that 
the eigenmode of the diffusion equation is a spin configuration exponentially decreasing in 
time, but its shape remaining unchanged. Taking the second order correction in w/lsp, one 
obtains 

i- M,k,m — i-M,k,m~^ TTa 72 l^iOaj 

- r° + 2(2 - M^)wWl 

for I A; — M/lsp ^/IsPi where r^o = Igp/vpl and Isp = Lso/2 have been used in the 
second equality. The modes with m = 0, and k = M/lsp will relax most slowly, since in 
these cases the first term r°^^^ in Eq. (15) disappears. The second term indicates that 
the spin relaxation time = I/Fm,*;,™ will be proportional to l/w"^. In the limit Lso ^ w, 



we thus 
systems. 



lave Ts much larger than the D'yakonov Perel' relaxation time r^o in 2D boundless 
2! Such a behavior will become more clear from the following simple consideration 
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in a ID system. Let us consider an electron at xi in a ID channel with an initial spin 
pointing to rii, as shown in Fig. 5(b). Under the Rashba SOI, the spin will rotate to n2, 
71-3, when this electron moves to X2, X3, .... It is easy to see that if each electron spin in 
an initial spin density distribution follows this (xj, n-j) relation, such a distribution will not 
change in time. Hence, its relaxation time is infinite. In a realistic 2D wire the relaxation 
time is finite at finite width. That is because electrons there can move in the y-direction. 
The polarization Pz{t) contributed from electrons moving along the channel is frozen, as in 
the ID case, while electrons moving along the the y-axis give rise to the relaxation of the 
helix distribution. 

Since the above expressions have been obtained under the assumption / ^ w ^ Lg^, it 
is interesting to extend the analysis beyond this limits. Within the SPI method we studied 
the relaxation of the helix mode in the range / ^ w, choosing Lgo = 12.5 /im, I = 0.5 fim, 
and vp = 0.37;um/ps. Given the initial helix spin mode \E'_|_i — \E'_i corresponding to spins 
oriented along the z-axis at x=0 in Fig. 5(a), the spin relaxation time Ts{w) calculated by the 
SPI method is represented by squares in the inset of Fig. 6. In the range of u; < 10 fim, this 
time increases dramatically and strongly deviates from the Tg of the corresponding uniform 
mode (red solid curve). Note that in order to display the divergent r^, we choose a large 
r^-axis scale in the inset of Fig. 6. At this scale the uniform mode (red curve) almost 
overlaps with the = axis. The relaxation time of the uniform mode will saturate at 
some value for w —>■ 0. It is ~ 28 ps for [Lso, I, w] = [2.19, 0.28, 1.4] fim in Fig. 1 and 
Ts ~ 441.1 ps for [Lso, I, w] = [12.5, 0.5, 0.4] /im. In contrast to these uniform modes, of 
the helix mode strongly increases for w (if w ^ /). This behavior is consistent with 
the theoretical result Eq. fll5ap . On the other hand, if w is as large as 20 /im, the relaxation 
time ~ 1927.4 ps of the helix mode is still much larger than ^ 241.0 ps corresponding to 
the uniform mode. This difference can be easily seen by magnifying the r<(-axis of the inset 
in Fig. 6. 

The main plot in Fig. 6 shows the l/r^ dependence on w"^. The squares show the helix 
relaxation rate corresponding to the data in the inset. The line with circles for the helix 
mode is obtained from the analytical expression Eq. fll5b|) . while the line with triangles is 
obtained from Eq. fll5al) . In the latter line, is simulated numerically, instead of using 
the analytical relation r^o = I'sp/vpl mentioned below Eq. fll5b|) . The lines with circles and 
triangles are valid only aXl <^w <^ Lso-, while that with squares is valid for all w. One can 
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clearly see that the curve calculated by the SPI method agrees very well with Eq. (15) at 
small w between wi = 2 fim (the smallest calculated width) and W2 = 10 fim. It is interesting 
to note that although formula Eq. (15) was derived under the condition / ^ w ^ Lgo, it 
seems to be valid in a wider range of w, because Wi is close to the ballistic regime crossover 
point at w ~ / = 0.5 fim and W2 is close to L^o = 12.5 /im. Finally, one should not expect that 
the linear trend in Fig. 6 can continue down to ^ because in this range the system will 
eventually reach the ballistic regime and Pz{t) will decay like a crossover function between 
the exponential and the Bessel functions. Similar to the discussion of the uniform initial 
spin configuration in Section III, it does not make sense to consider at extremely small 
w, because Pz(t) is no longer an exponential function. 



VI CONCLUSION 



The semiclassical path integral method has been applied to study the spin relaxation in 
thin 2D wires with the Rashba spin-orbit interaction. We considered the relaxation of a 
uniform spin polarization along the well as of the long-lived helix mode. In the 

former case we found a good agreement of calculated in the regime of large w {w ^ 20 /xm) 
with the well known bulk DP spin relaxation rate and with the experimental data from Ref. 
Q. At smaller w our numerical results deviate slightly from the experimental data. The 
nature of this distinction is not clear. We assume that the edge spin diffusion modes can 
contribute to the spin relaxation, so that the conditions for electron reflections from the wire 
lateral boundaries become important. Also, the Dresselhaus spin-orbit interaction can give 
rise to the observed dependence of Ts on the orientation of the wire axis in the xy-plane. At 
w ^ the relaxation time has a tendency to saturate at a value which is about twice of 
the bulk Tg, as predicted by the spin diffusion theory. Although both SPI and experimental 
data show similar saturation behavior, one must take into account that at w approaching 
the crossover w ^ If to the ballistic regime the evolution of the spin polarization can not 
be described by an exponential function. Hence, in this regime is not a representative 
parameter to describe the spin relaxation. We studied the evolution of the spin polarization 
in the ballistic regime and found that it is described by the Bessel function. The numerical 
SPI results fit well to this behavior. For the helix spin distribution, the linear dependence 
of 1/ts on w"^ predicted in the framework of the diffusion theory ^ coincides precisely with 
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that calculated by the SPI method. The SPI method also allowed to calculate the spin 
relaxation beyond the constraints of those analytic results. j4| 
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FIG. 1: (Color online) The spin relaxation times r^'s of [1,0,0] sample (circle) and [1,1,0] sample 
(star) versus channel width w are taken from the experiments in Ref. Q]. The calculated by the 
SPI method are extracted from the polarization curve Pzit), fitted by Eq. ^ with free parameters 
A and c (red solid curve) and with fixed parameters A = 1 and c = (black dash-dotted curve). 
The experimental Ts saturates at 11.5 /^m (pink dashed straight line) for large w, the same as the 
analytically estimated value for Lgo = 2.19 /um. The inset demonstrates three examples of Pz{t) 
for channel width w = 1.6 /im, 2.8 ^m, and 4.4 ^m. 
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FIG. 2: (Color online) In the inset, a spin configuration in a channel of the length Svr fim relaxes 
to zero, depicted at three different times. These configurations are spatially uniform up to the 
ripples at two ends caused by boundary effect. Recording the polarization Pz{t) at the middle 
point of the channel gives the relaxation curve (green solid thick) in the main plot. This curve 
can be fitted by the exponential function in Eq. ([6|) with [r^, A, c] = [7.127, 1.274, —0.055] (blue 
dashed curve) and [rs,A, c] = [4.323, 1, 0] (red dotted curve). For t close to zero, most electrons 
have not been refiected by impurities or boundaries. In this range Pz{t) does not behave as an 
exponential function. Later, after most of the electrons and their spins have been randomized by 
impurities or boundaries, Pz{t) became more exponential-like. The physical parameters used are 
[w, Lso,vf,1] = [0.1 /im, 2 /xm, 0.37;um/ps, 0.3 /im] with 6 x 10^ electrons. 
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FIG. 3: (Color online) In a ID channel without impurities, a spin polarization Pz{t) behaves like 
a sinusoidal function Eq. ([8]) (blue dash-dotted curve). In an infinitely thin channel without 
impurities, Pz{t) behaves like a Bessel function Eq. ^ (smooth red dashed curve), which agrees 
with the numerically obtained Pz{t) simulated by 5 x 10*^ electrons (rugged green solid curve). The 
physical parameters are [w, Lso,vf,1] = [0.1 ;um, 2 /xm, 0.37/im/ps, 10^ /im]. 
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FIG. 4: (Color online) Two kinds of trajectories with the same length I = vpt in a thin channel: 
the straight trajectory pTpo and the zigzag trajectory p2Po- The latter has the same length as the 
straight line P2P3, since it can be obtained through a multiple mirror reflection of P2P3 with respect 
to the horizontal dotted lines. The symbol 9{^) denotes a certain outgoing angle of an electron 
located at ^, as explained in the text. The trajectory P2PQ is magnified in the inset. 
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FIG. 5: (Color online) (a) Long-lived spin eigenmodes: ^'o denotes the spin mode with all spins 
aligned in y-direction. — ^'^i represents the helix spin modes with spin rotating on the xz- 
plane. (b) An electron in the channel moves a distance ttLso from the left end xi to the right end 
xg. Due to Rashba SOI, the spin of the electron precesses from ni to ng and completes a phase 
period of 27r. (c) A schematic plot of the decay of the helix mode at 1 (dotted), 35 (dash), and 170 
(solid) time units. "^±1 in Fig. 5(a) has been replaced by — ^'-i. 
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FIG. 6: (Color online) The relaxation time of the uniform mode (red solid) and the helix mode 
(black square) versus the channel width w in the inset. For small w, the square curve of the helix 
mode is redrawn as l/r^ against w'^ in the main plot, to compare with the r^'s derived from the 
analytical formula (|15bp (red circle) and the semi-analytical formula ()15ap (green triangular). 
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